Electron— electron Interactions in Disordered Metals: Keldysh Formalism 
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We develop a field theory formalism for the disordered interacting electron liquid in the dynamical 
Keldysh formulation. This formalism is an alternative to the previously used replica technique. In 
addition it naturally allows for the treatment of non-equilibrium effects. Employing the gauge in- 
variance of the theory and carefully choosing the saddle point in the Q-matrix manifold, we separate 
purely phase effects of the fluctuating potential from the ones that change quasi-particle dynamics. 
As a result, the cancellation of super-divergent diagrams (double logarithms in d = 2) is automati- 
cally build in the formalism. As a byproduct we derive a non-perturbative expression for the single 
particle density of states. The remaining low-energy cr-model describes the quantum fluctuations 
of the electron distribution function. Its saddle point equation appears to be the quantum kinetic 
equation with an appropriate collision integral along with coUisionless terms. Altshuler-Aronov 
corrections to conductivity are shown to arise from the one-loop quantum fluctuation effects. 

PACS numbers: 71.10.Ca, 71.23.An, 72.10.Bg. 



I. INTRODUCTION. 



The physics of weakly disordered interacting electron 
systems at low temperatures has been a subject of con- 
siderable theoretical and experimental interest over the 
past years (for review see Refs. |^J^ ). Although sig- 
nificant progress has been made in this direction, many 
thermodynamic and transport properties of such systems 
are not completely understood and continue to stimulate 
both experimental and theoretical research. The latest 
revival of interest to the problem was prompted by the ex- 
perimental discovery of a possible metal-insulator tran- 
sition originally in clean Si MOSFETS Q and later in 
and p-type GaAs 

The low temperature behavior of the conductivity of 
a metal is mainly determined by the quantum (weak 
localization) jsl^ and the interaction corrections to 
the classical Drude result. These corrections are espe- 
cially strong in low dimensional {d < 2) systems. In 
two-dimensions, for example, both the lowest order weak 
localization correction |^ and the lowest order interac- 
tion correction Q diverge logarithmically at low temper- 
atures. The ultimate faith of the low temperature phase 
is determined by the interplay between them. 

According to the scaling theory of localization Q , in 
the absence of electron-electron (e-e) interactions (and 
with no spin-orbit scattering), the quantum corrections 
lead to localization of all single particle states in dimen- 
sions d < 2 and thus to insulating behavior for arbitrarily 
weak disorder (weak localization) . Wegner |^ proposed a 
replicated cr-model to study this problem. With the cou- 
pling constant corresponding to the dimensionless con- 
ductance g, this cr-model provided justification for the 
one-parameter scaling theory of localization |^]. Later, 
Efetov [ p^ introduced the supersymmetric version of the 
cr-model which obviated the need to take the tricky pll] 



zero replica number limit. 

Finkel'stein jl^ developed a replicated cr-model ap- 
proach for interacting disordered electron systems, which 
was farther developed in |l3Ul^ . He demonstrated fur- 
ther its renormalizability in the one-loop approxima- 
tion and obtained the one-loop renormalization group 
flow equations. From these equations it followed that 
the weak coupling fixed point corresponding to non- 
interacting metal is unstable. The need for introduc- 
ing replicas in Finkel'stein's approach follows from the 
fact that the ensemble averaged observables are obtained 
as derivatives of the averaged logarithm of the partition 
function. The formalism in Ref. JT^-p^ utilizes Mat- 
subara representation and is therefore restricted to the 
equilibrium situation. 

Later it was suggested that the Keldysh type field 
theory, originally developed for the treatment of non- 
equilibrium systems [ [l5[ may be an alternative to the 
replica technique [p6[-p^. The point is that the use of 
the Keldysh closed contour in the time direction, leads 
to an automatically normalized (disorder independent !) 
partition function. This circumvents the need to intro- 
duce replicas. A similar situation exists in the theory of 
spin glasses, where in addition to the replica approach 
p9[ the Martin-Siggia-Rose formalism , analogous to 
the Keldysh approach [p0|-p2[, has been used. This for- 
malism provided insight complementary to that gained 
from the replica approach. Horbach and Schon |^ de- 
veloped a cr-model for non-interacting electrons in the 
Keldysh formalism. Although our treatment differs from 
their in many important details, we have benefited much 
from their work. 

Here we apply the Keldysh formalism to disordered 
interacting systems. We restrict ourselves to the consid- 
eration of spinless electrons in the presence of a weak 
magnetic field (unitary ensemble) and leave the con- 
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siderations of the spin and Cooper channels for future 
work. Another important distinction of the present the- 
ory from the previous ones ]T^-p^ is the different choice 
of a saddle point of the functional integral on the Q- 
matrix manifold. The saddle point in our formalism ex- 
plicitly depends on a fluctuating potential in the sys- 
tem (the Hubbard-Stratonovich field, which decouples 
the electron-electron interaction). This choice of the sad- 
dle point allows us to separate pure phase effects of the 
fluctuating potential and to present the first, in our opin- 
ion, clear derivation of the tunneling density of states 
(DOS) in a metal film obtained earlier by Finkel'stein 
[g2| and Levitov and Shytov |Q by different means. An- 
other advantage of this choice of the saddle point is that 
the perturbative expressions for gauge invariant quan- 
tities contain only single logarithms of temperature or 
frequency (in d = 2). The diagrams containing dou- 
ble logarithms which appear in the standard diagram- 
matic expansion |25| or in Finkel'stein's formalism fl^ ] 
and cancel each other for any gauge invariant quantity 
do not appear in our formulation at all. This signifi- 
cantly reduces the number of diagrams in each order of 
the perturbation theory. We then obtain a low energy 
theory in the form of a cr-model. The advantage of the 
Keldysh formulation is that it allows for a clear physical 
interpretation of the effective degrees of freedom. They 
turn out to be the quantum fluctuations of the electron 
distribution function. The saddle-point equation on the 
massless manifold is just the quantum kinetic equation 
with an appropriate collision integral. The one-loop fluc- 
tuations on top of this saddle point lead to corrections to 
various observables and in the case of conductivity can be 
identified with the Altshuler-Aronov corrections [[Zp^]. 

The paper is organized as follows: In section |l| we 
present the functional integral representation of the 
Keldysh partition function for disordered interacting 
electron systems. Section III is devoted to the choice 



of an interaction-dependent saddle point and the deriva- 
tion of an effective tr-model as the massless fluctuations 
around this saddle point. We discuss some applications 
of the theory, like the derivation of the non-perturbative 
expression for the single-particle Green function, in sec- 
tion IV. Quantum fluctuations and Altshuler-Aronov 



corrections to the conductivity are the subject of section 
In section VI we derive the quantum kinetic equation 
as the saddle point equation on the massless manifold 
and discuss the corresponding collision relaxation times 
along with the collisionless terms. Finally, in section VII 
we briefly discuss the obtained results and the future per- 
spectives. 



II. FUNCTIONAL INTEGRAL FORMULATION 



A. Keldysh formalism 

Consider a unitary evolution of a system along a closed 
contour C in the time direction which consists of the prop- 
agation from t ~ ~oo io t = -\-oo and then back from 
t — -foo to t = — oo. All external time-dependent fields 
are assumed to be exactly the same during the forward 
and backward evolution processes. As a result, at the end 
of such evolution the system must find itself precisely in 
the original state. We thus conclude that the evolution 
operator 

Uc = l. (1) 
Let us consider next the partition function defined as 



Z = Tr{po?/c}/Tr{po} = 1 ■ 



(2) 



where po is a density matrix of the system at the ini- 
tial time, t = — oo, before the interactions and disorder 
are adiabatically switched on. A more informative ob- 
ject is the generating functional, which is obtained by 
introducing source fields. It is clear that to have a gener- 
ating functional not identically equal to unity, the source 
fields should have a different behavior on the forward 
and the backward parts of the contour. To shorten the 
subsequent expressions we shall operate with the parti- 
tion function, Eq. (|^), and will introduce the generating 
functional in section |1V C , 



- oo +CO 
FIG. Schematic representation of the discretization of the 
time contour C. The dots on the upper and the lower branches 
of the contour denote the discretized time points. 

The next step is to divide the C contour into 2N + 1 
time steps, such as ti = t2N+i = —oo and tAr+i = -foo 
as shown in Fig.|l|. Following the standard route we 
obtain the coherent state functional integral, by intro- 
ducing a resolution of unity at each time step. Taking 
the N oo limit we obtain for the partition function 



Z = N X>VV'exp{iS'[?/',?M}, 



(3) 



where J\f is disorder independent normalization constant 
p8[ and the fermionic action is given by 

S[^, V] = y dt\ j dvii[G^^ ~ Ud^s{v)\i^ - (4) 
c ^ 

^ ^ ^ drdr'^(r)^(r')Vb(r - r')V'(r')V'(r) 

Here the inverse bare Green function is a shorthand no- 
tation for 
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Go ^ = z 



(5) 



d_ Yl 

dt ^ 2m ' 

where the time derivative is taken along the contour C. 
The notation is somewhat symbohc: while inverting 
this operator it is necessary to invert its discretized ver- 
sion first, and only then take the limit N ^ oo |^7|] . 

We then divide the fermionic field ip{r,t) into the two 
components '>pi{r,t) and "02 ("''j^) which reside on the for- 
ward and the backward parts of the time contour respec- 
tively. Since the interaction part of the action is strictly 
local in time it may be rewritten as Smtii^i] — Sint['>p2\ 
(the minus sign comes from the opposite direction of the 
time integral on the backward part of the contour) 

S^nt m = ^Jdt dr dr% (r)^, (rO^o (r - r') (r') (r) , 

(6) 

here i = 1,2 and Vb(r — r') is a bare interaction po- 
tential. We now introduce two independent auxiliary 
bosonic fields (f>i(2) (r, t) to decouple the two interaction 
terms by the Hubbard-Stratonovich transformation. As 
a result one obtains for the partition function 

$] =Tr {W[Go 1 - Ud,s<y3 + (^aT"]*} • (7b) 

Here we have introduced the following vector notations 
for the fermionic doublet the bosonic doublet $, and 
two vertex matrices 7" 



* = 



V'2 



$ = 



(8a) 



7 



1 




7 





-1 



(8b) 



The inverse matrix Green function stands for 



Gq 



-1 



"dt ^ 2m 







; J. _ Zi 

' dt 2m 



(9) 



The trace operation in Eq. (7b) and henceforth is under- 
stood to be performed over the 2x2 structure as well as 
over the time and space variables. 



B. Disorder averaging 

The great advantage of the Keldysh technique is that 
the normalization constant, J\f in Eq. ( [7a| ) does not de- 
pend on the realization of the disorder potential. Thus, 
the disorder averaging can be performed directly, with- 
out the need to resort to the replica trick. Hereafter we 
employ the simplest model of the Gaussian, (5-correlated 
disorder 



y I?f/d»s • ■ -exp |-7rz/r J drUi,^{r)^ 



(10) 



where the disorder strength is characterized by the elastic 
mean free time, r; v is the bare density of states at the 
Fermi energy. Next, we perform the Gaussian integration 
over Udis in Eq- ( [7a| ) and decouple the arising (non-local 
in time) quartic interaction by means of the Hubbard- 
Stratonovich transformation. Doing so, we obtain 



jTr{*c/<iisO-3*} 



exp ■ 



J 2?(5 exp I - y drdtdt' 



dr 



dt'J'(r,t)(73*(r,<) 



4r 2t 



(11a) 



(lib) 



Here we have introduced the Hubbard-Stratonovich field 
Q which is a matrix with indices both in the Keldysh 2x2 
space and in the time space. To ensure the convergence 
of the integral in Eq. ( [llb| ) the Q-matiix is chosen to 
be Hermitian. After these transformations the fermionic 
functional integral in Eq. ( |7a| ) can be formally performed, 
leading to 



det 



Gn 



2t 



+ (/'a7° 



(12) 



As a result, the disorder averaged partition function takes 
the form 



(Z) 



(13a) 



iS[Q,'l>] = TrQ2 + Trln 

4t 



Gn 



2t 



*a7 



(13b) 



As before the trace operation is understood to be per- 
formed over the Keldysh and the time indices as well as 
over the coordinate space; the unessential normalization 
constant is omitted. 



C. Keldysh rotation 

In the notations introduced in Eqs. (H) and (|^) the elec- 
tron Green functions G are matrices in the 2x2 Keldysh 
space. Their components are not independent and satisfy 
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certain general identities ||l^,g9[. This interdependence 
becomes most transparent if one introduces the rotated 
Green functions G denoted by the absence of the hat and 
defined as 



G EE LcrsGL'' , 
where the unitary matrix L is given by 



(14) 



(15) 



As follows from the definition of the Keldysh Green func- 
tion [|l^,^ , the rotated Green function has the following 
structure 

where G^^'^Ht^t') vanish for t < t'{t > t'). To pass to 
the rotated representation we introduce a new Hubbard- 
Stratonovich field Q which is related to the old one, Q, 
by the following unitary transformation 



= LQU 



(17) 



We also introduce the rotated bare inverse Green func- 
tion, Gq^ , expressed through Gq ^ of Eq. (^ in a manner 
consistent with Eq. (|lj) 



Gq ^ = LGq ^(TsL^ 



. d 



(18) 



It is also convenient to perform a linear transformation 
of the bosonic doublet, $, by introducing the symmetric 
and the antisymmetric combinations of the fields residing 
on the upper and the lower branches of the contour C 



h = i(0i + 02); 

h = h{(t>i - 02) • 



(19) 



Then the rotated vertex matrices for these new fields are 



7 



1(2) 



L(7^ ± 7^)ct3L^ with the following explicit form 



7 



1 
1 



7 



1 

1 



(20) 



Any classical external field takes on identical values on 
the two branches of the contour and, hence, in the ro- 
tated basis has only the first symmetric component. The 
non-zero antisymmetric component may appear only as 
a virtual fluctuating field. Below we shall sometimes refer 
to the first and second components of the bosonic fields 
as classical and quantum ones correspondingly. Since the 
presence of an external classical field does not change the 
basic fact that Z — 1, any auxiliary source field should 
have a non-vanishing quantum component to generate an 
observable. We shall return to this observation in section 



IV C 



Utilizing the cyclic invariance of the trace operation we 
obtain the following expression for the partition function, 
Eq. (|l3|), through the new variables $ and Q 

{Z) = I P$e'Tr{$-y-Vi*} J 'pQ^^S[QM ^ (21a) 

i5[Q,$] =-^TrQ2 + Trln[Goi + :^Q + 0a7H. (21b) 
4r L zr J 



III. NON-LINEAR cr-MODEL 

A. Saddle point equation 

We shall look now for a saddle poi nt of the functional 
integral over the Q- matrix in Eq. (21a). The aim is 
to find a stationary solution for a given realization of 
the slowly varying in space and time fluctuating b osoni c 
flelds, Calculating a variation of the action, Eq. ( |2Hj ), 
over the Q-matrix, one obtains the following equation for 
the saddle-point matrix, Q = 



) (r) = — 



Gq ^ 



+ :rQ + 0a7° 

ZT 



(22) 



We are unable to solve this equation exactly, therefore 
our goal will be to find its approximate solution, which 
is as close as possible to t he t rue stationary point of the 
functional integral in Eq. ( 21a ). To execute this program, 
we first consider the case where $ = 0. It is easy to check 
that in this case 



A*_i, =Q[<i> = 0] = — VG(p,t-i'), 



(23) 



where the impurity averaged Keldysh Green function is 
(cf. Eq. ^) 



G{p,e) 



1 

-1 



G^ip,e) G^ip,e) 
G^{p,e) 

G«(p,e) 
G'^(p,e 



(24) 



1 

-1 



with 



G^(^)(p,e) = (e-ep±*/(2r))-i 
G^(p,e) = G^F-FG^. 



(25a) 
(25b) 



The function F defined by Eq. (25b) can be expressed 
through the single particle distribution function, n(e), as 
F{e) = 1 — 2n(e). In equilibrium at temperature T it is 
given by 



tanh — - 
2T 



(26) 



Substituting Eqs. (|2|), (|5|) into Eq. (|3|) and perform- 
ing the momentum summation, one obtains for the non- 
interacting ($ = 0) saddle point 
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A, 



af 





A, 







(27) 



We have introduced here the retarded and the advanced 
unities, if which should be understood as Fourier 
transforms of infinitesimaUy shifted 5-functions. This 
particular form of the Green function is a result of the 
approximation that the single-particle DOS is indepen- 
dent of the energy, e. In reality it does depend on e, and 
the retarded (advanced) components of A(e) are analytic 
functions of energy in the upper (lower) half plane which 
do depend on energy on the scale of order of the Fermi en- 
ergy, tp- Therefore the infinitesimaUy shifted (5- functions 
in Eq. ( p7| ) should be understood as St^i^ = f±{t)0{±t), 
where Q{t) is the Heavyside function, and f±{t) are func- 
tions that are highly peaked for \t\ ^ ej,^ and satisfy the 
normalization condition °° dtf± {t) — 1 . This high- 
energy regularization is important to remember in calcu- 
lations to avoid spurious unphysical constants. In partic- 
ular, for obvious reasons 



R 

t-t' 

A 

t-t 



'MP,, = 



0, 



(28) 



where M'^^f'^ is an arbitrary retarded (advanced) matrix 
in the time space. 

Substituting Eq. (|2^) into Eq. (|2|) with $ = 0, it is 
easy to see that Q = A solves the non-interacting saddle- 
point equation for any function F^. (The simplest way 
to check it is to use the decomposition Eq. (^).) This is 
natural, since any distribution function is allowed for the 
non-interacting electron gas. We shall see below how the 
interaction effects drive the system towards the equilib- 
rium distribution, Eq. (p6|). 

Let us now include a finite 4>a {r, t) into Eq. ( |2^ ) . To 
this end we notice that this equation can be still solved 
exactly for the particular case of spacially uniform real- 
izations of the boson field, (pa — 4>a{i)- This is obvious 
since such a field may be gauged out resulting in 



(29) 



The validity of this solution can be verified by acting with 
the operator [Gq ^ + i/(2T)(5 + (f)a{t)^°'] on both sides of 
Eq. ( p^ ) and utilizing the fact that Kt-t' solves Eq. ( |2^ ) 
with $ = 0. We also rely on the commutativity of the 
vertex matrices, [7^,7^] = 0, in writing the solution in 
the form of Eq. (|29|). 

We now consider the case where the bosonic fields (pa 
are slowly (compared to the mean free path I) varying 
in space. In analogy with Eq . (|29| ) we shall look for an 
approximate solution of Eq. (|22|) in the form of a local 
(in time and space) gauge transformation of A 



At-t'e 



(30) 



where ka = is a certain linear functional of the 

fields $, whose specific form is to be determined to sat- 
isfy Eq. (HI) in the best possible way. 

To proceed we introduce a new Hubbard-Stratonovich 
field Q which is related to the old one, Q, by the gauge 
transformation 



-ikct (r,i')7° 



(31) 



Substituting this definition into the action, Eq. (20), and 
using the invariance of the trace under a cyclic permuta- 
tion of operators, we can rewrite the action as 



*5[Q,a>]--— TrQ2 



(32) 



Tr In 



(Vfcg) 

2m 



where we have introduced the notation 

C(r, t) EE (0„ - dtka - ^pVka)!'' , 



(33) 



with the Fermi velocity, 'Wp = — zV^/m. To find the ap- 
proximate saddle point of the form Eq. ( pO| ) we substitute 
Q = K + SQ into Eq. ( ^2| ) and require terms linear in SQ 
to vanish. In doing so we neglect the diamagnetic term, 
(Vfca)^/2m, since it is quadratic in ka (and hence in $) 
and is also smaller than C in the parameter q/pp <^ 1, 
where q is the characteristic momentum scale of variation 
of <&. As a result we obtain the following equation 

- TTi'At'.t + i[G- GCG + GCGCG - . . .],, , (r, r) = . 

(34) 

The first two terms in this expression cancel, according 
to Eq. (p3|). The freedom of choosing the K[^] functional 
is not sufficient to cancel all the terms in this expansion. 
We thus concentrate on the term which is linear in <I> and 
K: 

Y,G{p+,e+)C{q,Lu)G{p-,e-) = niyr x (35) 
p 

[((/.a+zc^fca)(7"-A+7"A_)-i?g2^„(A+7"-7"A_)] , 

where p± — p ± g/2, e± = e ± uj/2 and A± = A^^ . To 
derive Eq. ( ^5| ) one may employ the following useful rep- 
resentation of the Keldysh Green function 



G(p,e)^[Goi + -A,]-i 

= iG«(p, e)(ao + A,) + ^G^{p, 6)(ao - A,) 



(36) 



Only J^pG^G^ and "^G^vpG"^ contribute to 
Eq. (p|). Multiplying Eq. IM) by A^^^ from the left one 
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obtains the following matrix condition for the vanishing 
of the linear term in Eq. ( ^4| ) 

Dq^k^ (A+7"A_ -7") + (0„ + zu;fc„) (A+7" -7" A_ ) = 0. 

(37) 

To cancel (1, 1) , (2, 2) and (2, 1) components of the ma- 
trix on the l.h.s. of this equation the functional K should 
satisfy 



(38) 



Provided this equality is obeyed, the condition to cancel 
the Keldysh (2, 1) component on the l.h.s. of Eq. (|3^) is 



1 — F F 

(Dq^ - iuj)ki -(t>i^ -2Dq^k2 



F, 



F,_ 



(39) 



This equation can not be in general satisfied, since its 
r.h.s. contains explicit dependence on e, whereas the 
l.h.s. is supposed to be e-independent. This happens, 
because the trial saddle point solution, Eq. (30), is too 
restrictive. In particular, we demanded that it may be 
obtained from A by the rotation which is local in time. 
As a result, the functional K depends only on differen- 
tial energy uj and not on "center of mass" energy e. This 
restriction is in an apparent contradiction with Eq. (39). 
There is, however, an important particular case when 
Eq. (^9|) may be solved. This is the case of thermal 
equilibrium, where the fermionic distribution function is 
given by Eq. (Eq). As a result 



F!l - F!^_ " 2r = ' 



(40) 



where stands for the equilibrium bosonic distribu- 
tion function. Thus, the matrix equation (^) for the 
functional K[<i'] can be resolved in equilibrium. The re- 
sult may be written in a short form as 



V-\q,uj)K{q,uo) ^liZ^^iq.^) ■ 



(41) 



We have introduced here the following bosonic matrix 
propagators in the 2x2 Keldysh space 



V{q,uo) 



V^{q,Lu) V^iq,co) 
V^{q,cj) 



with 



V^'^^\q,cu)^{Dq'T^^ 



V^{q,L0)^B^{V'\q,u)-V''{q,u)): 



and 



'iujV{q — 0, uj) 



25. 
-l^, 



(42) 



(43a) 
(43b) 



(44) 



The superscript "eg" denoting equilibrium in the bosonic 
distribution has been omitted for brevity. 



Eqs. ( |30| ) and (41) complete the task of finding the 
approximate saddle point, Q = Q[^], for any given real- 
ization of fields <I>. On this solution we are able to cancel 
only the term linear in $ in the expansion Eq. ( ^4[ ) . This 
guarantees only that the terms like ^5Q will not appear 
in the expansion of the action around the saddle point 
given by Eqs. (|^), (|l|). Terms like ^^5Q may (and 
will) arise in such expansion. We shall see later, that 
it is precisely these terms that are responsible for the 
divergent Altshuler-Aronov corrections to conductivity 
p6[ . The ability to avoid '^5Q terms is, strictly speak- 
ing, limited only to the thermal equilibrium. For an out 
of equilibrium situation such terms reappear and require 
some care (see Section VI). 



The influence of the external potential <& on the elec- 
tron dynamics (and hence on the Green function) is two- 
fold i) It changes the particle trajectory and ii) It 
changes the phase of the electron wave function. The first 
effect is proportional to the electric field E = — V$ and 
is small for the long wave length spatial configurations of 
<I>. The second effect, however, requires no actual elec- 
tric fields. It is proportional to $ itself, rather than V$, 
and is akin to the Aharonov-Bohm effect. It changes the 
phase but not the amplitude of the wave function and can 
be taken into account in the Eikonal approximation. The 
second effect exceeds the first one for the long wave length 
fluctuations of the potential, therefore, it is especially im- 
portant in the presence of the long range Coulomb inter- 
actions. The approximation to the saddle point, Eq. (|3^), 
is similar to the Eikonal approximation. It is designed to 
account for the phase effect of the slow fluctuations of 
the potential $. Note that the phase K enters the sad- 
dle point equation only through its total time derivative 
along the trajectory of a particle, d/dt = dt + Vi?V (cf 
Eq. (^)). If we demand that C vanish we obtain the 
standard Eikonal equation for the action K of the 
particle moving with a given velocity v in an external 
field $. Unfortunately, the ansatz Eq. (|3C| ) is too re- 
strictive to nullify C for particles of every velocity v. 
Eventually all the particles in the Fermi sea interfere to 
produce the Green function Q^^X^)- Equation ( ^o|) ap- 
proximately accounts for the phase interference between 
particles moving along different trajectories. Since the 
particle dynamics is diffusive this leads to the diffusive 
relation (^l|)-(^3|) between the external potential $ and 
the phase K. As will be clear below the choice of the 
saddle point in the form Eq. (^0|) considerably simplifies 
the subsequent calculations. In particular, it eliminates 
completely the family of super-divergent diagrams which 
cancel in the traditional treatment after sometimes 
tedious calculations. 
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B. Effective action 

To formulate an effective low energy theory in terms 
of the fluctuating fields Q and <i> we need to examine 
the fluctuations around the saddle point. The fluctua- 
tions of Q fall into two general classes [p| jlO|Jl^ : i) mas- 
sive, with the mass oc 1/r and ii) massless, those on 
which the action depends only very weakly. The fluc- 
tuations along the massive modes can be integrated out 
in the Gaussian approximation and lead to insignificant 
renormalization of various parameters in the action. The 
massless, or Goldstone, modes describe diffusive motion 
of the electrons. The fluctuations of the Q-matrix along 
these massless modes are not small and should be treated 
carefully. The Goldstone modes can be parameterized by 
the Q-matrices satisfying a certain nonlinear constraint 

MB- 

To identify the relevant Goldstone modes consider the 
flrst term in Eq. (32). The saddle point is given by 
^) and (|4^) satisfies 

<3^ = ( ) . m 



li 



Eqs. 



and the first term in Eq. ( |32D vanishes. The fiuctuations 
of Q which do not satisfy Eq. (^) are massive. The 
massless modes are generated by rotations of the saddle 
point and can be parameterized as [^,^,|2| 



T- 



(46) 



The parameterization of the rotation matrices T must 
ensure the convergence of the functional integration over 
the matrices given by Eq. (^6|). Below we only assume 
that such a parameterization exists, whereas the concrete 
form of T is not important for what follows. 

One way of parameterizing the rotations is to write 
T = exp{M^/2}, where, without loss of generality, WA = 
—AW. Expanding Eq. ( |3^ ) to the second order in W and 
neglecting for a moment the term arising due to e-e inter- 
actions it is easy to establish that in the diffusive regime 
the relevant fluctuations must satisfy the condition 



W,,,,^0, only if |e|,|e'| < l/r. 



(47) 



Namely, all effective degrees of freedom are concentrated 
in the narrow energy strip of the width 1/t <C cf 
near the Fermi energy. Therefore the matrices T dif- 
fer from unity only in the narrow region of energies de- 
flned by Eq. (^^. For this reason the gauge transfor- 
mation Ut^t'ir) = ex.p{—ika{t, r)ja}S{t — t') in Eq. (p^) 



can not be incorporated into a redeflnition of T and 
should be carried out explicitly. Indeed, being diagonal 
in time indices, the matrix hlt,t' spreads over the entire 
energy space and, thus, can not be reduced to a distur- 
bance, which is close to the Fermi shell. Physically, this 
describes the fact the low-wavenumber scalar potential 
<i>(g, t) shifts the entire electronic band and not only the 
energy strip given by Eq. (^^. It is essential to follow 
the variations of the electron spectrum all the way down 
to the bottom of the band to respect the charge neutral- 
ity imposed by the Coulomb interactions. As we shall see 
below, once the phase factors in Eq. ( p6| ) have been taken 
into account, the residual interactions may be regarded 
as being short-range without loss of generality. 

Substituting Eq. (^) into Eq. (|3^), and retaining only 
the universal (r - independent) terms in the expansion 
of the logarithm, we obtain for the Q action 



wTt{{^ -iujKfaii^ + iLjK)} ~ (48) 
DTridrQY + 4zTr{(e + (0„ + iwfc„)7")g}l , 



where we have introduced the long derivative 



(49) 



A few comments are in order regarding Eq. (|4^). First, 
it is restricted to Q which satisfy Eq. (^5|). The last 
two terms, containing Q, conventionally originate from 
^pVpG^vpG^ and ^^G^^"^^ combinations in the ex- 
pansion of the logarithm. On the other hand, the flrst 
term in the r.h.s. of Eq. ( ^ ) originates from J^p G^G^ 
and G^G^ combinations. These terms should be 
retained since, as was mentioned above, the matrix 
0a (e — e')7" is not restricted to the l/r shell near the 
Fermi energy. To derive this term we employed the fact 
that for any physical fermionic distribution function 



*±oo 



±1. 



(50) 



Finally, the terms like v^G^vpG^, although non- 
vanishing, cancel against the diamagnetic term. 

Employing the explicit form of the long derivative, 
Eq. (Eof), and the relation between the K and $ flelds. 



Eq. (41), one finally obtains for the the partition function 



{Z)= jv<^ exp{iTr{$^y-i$}} jvQ exp {i^oiQ] + VX] + ^52[g, V/f]} , 

where S';, ^ = 0, 1, 2 contain V K in the I — th power and are given by 



iSa[Q] =-— i^Tr{VQ}' + 4i Tr{eQ} 



i^Tr{Vfta7"gVQ} + Tr{(0„ + ia;fc„)7"Q} ; 



(51) 

(52a) 
(52b) 
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S2[Q, Vif] = -— Tr{VA:„7"QVA:;37^Q} - Tr{Vfc„7"AVA:/37'^A} 



(52c) 



The effective interaction matrix V is nothing but the 
screened interaction in the RPA approximation 



V{q,Lo) ^ {Vo-\q)ai + Po{q,io)) ' , 



(53) 



where PQ{q,uj) is the bare density-density correlator. It 
has a typical form of a bosonic correlator in the Keldysh 
space 



with 



r,R{A), ^ Dq^ 



(54) 



(55a) 
(55b) 



To derive Eqs. (pl|)-(|55|) we had to add and subtract the 
term Tr{VfcQ,7"AVfc/37''A} and employed the equation 



+ C30 

Jde Tr{7"7'3 - 7"A,^7''A,_ } = Acu {Uz') 



(56) 



Here e± = e±uj/2 and matrices A and 11 are defined by 
the Eqs. ( |27| ) and (^J) correspondingly. Eq. ( p6|) based 
on the following relations between bosonic and fermionic 
distribution functions 



Jde [F,^ - F,_) ^ 2u; ■ 

-OO 

+ 00 

Jde{l-F,^F,_) = 2cjB^ 



(57) 
(58) 



The last equation is obviously satisfied in the thermal 
equilibrium. For a non-equilibrium situation it should 
be considered as a definition of B^^. 

Eqs. (|5l])-(|5|) together with Eq. (|l]) constitute an 
effective non-linear cr-model for interacting disordered 
electron gas. The model consists of two interacting 
fields: matrix field Q, obeying the non-linear constraint 
Eq. (^), and the bosonic vector field <i> (or equivalently 
K). As will be apparent later, the Q-field describes 
fluctuations of the quasi-particle distribution function, 
whereas $ (or K) represents propagation of electromag- 
netic fields through the media. The following sections are 
devoted to the analysis of this model and calculation of 
various physical quantities on the basis of the model. 



IV. APPLICATIONS OF THE FORMALISM 

A. Single particle Green function 

In this section we shall show how the developed formal- 
ism can be used for the calculation of the average single- 
particle Green function at coinciding spatial points. This 
quantity is defined as 



(59) 



where ((. . .)) denotes both the quantum and the disorder 
averaging. It is convenient to apply the Keldysh rotation, 
Eq. (ift), and define 



g{t-t')^L<j3g{t^t')L^ 



(60) 



Such Green function arises e.g. in calculations of the tun- 
neling DOS, or shot noise power. To evaluate it one may 
introduce a source term in Eq. (jj), directly coupled to a 
bilinear combination of the fermion operators. Following 
the same algebra as above one finds that the source field 
enters into the logarithm in Eq. (2 lb). Differentiating 
finally with respect to the source and putting it to zero 
one obtains for the Green function 



(61) 



VQ e 



iSlQ,'S>] 



Zt 



We shall evaluate the integral over the Q matrix by the 
saddle point approximation, neglecting both the mas- 
sive and the massless fluctuations around the stationary 
point. Then, according to Eq. (p^, the pre-exponential 
factor is simply —iniyQ^ At the saddle point Q is given 
by Eq. (|30|). Transforming the action, S[Q, <&], in the way 
it was done in section IIIB, one obtains for the Green 



function in the saddle point approximation 
g = -»7r^'|l?<I)e'Tr{*-y-i*}g.fc„(th°A,_^,g 



(62) 



Since K is the linear functional of given by Eq. (|4l|), 
the remaining functional integral is Gaussian. Employing 
Eqs. (|l]) and (|5|)-(||), one obtains for the correlator of 
the K fields (averaged over the fiuctuations of $) 



{ka{q.,uj)ki3{-q,-uj))^ = -Vap{q,uj) ; 



(63a) 



V(q,L^)=I?(g,a;)njV(g,co)(n:i) V^{-q,-oj). (63b) 

The Keldysh matrix V has the familiar structure of a 
bosonic propagator 
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V{q,uj) 

with 



V^iq,uj) 



(64) 



-1 



1 



, (65a) 
(65b) 



V^(g,a;)=B.(V«(g,c.)-V^(g,t^)). 

One may recognize that this propagator precisely corre- 
sponds to the screened Coulomb interaction line dressed 
by two diffusons at the vertices. Thus, the role of the 
if-field is to take into account automatically both the 
RPA-screened interactions and its vertex renormalization 
by diffusons. 

To calculate the functional integral, Eq. (^), we write 
the phase factors as 



g±ifco7° _ _ |g±i(fci+fc2) g±i(fci -/C2) 



(66) 



±i(fci+A:2) 



and perform the Gaussian integration according to 
Eq. (|63|). The result may be conveniently expressed in 
the following form 



Git) = -^7^^. ^ {j'^Atl") B^p{t) . (67) 

Q/3=l 

where the fictitious propagator B has the standard 
bosonic structure (as e.g. Eq. (H^)) with 



(*) 



B'^it)^ iei(v'^(*)-v'^(o)) |^e*(v^(t)-v-*(*)) +e-*(^''(*)^'^'*(*») 



(68a) 
(68b) 



The {KK'^) propagator, V, defined by Eqs. (|6|),(||) is 
taken at coinciding spatial points 



Vit) 



(69) 



The electron Green function must satisfy several im- 
portant requirements: the tunneling DOS, i^(e), which is 
defined as 



Ke)-;^(^«(e)-t;^(e)) 



(70) 



must be a positive definite quantity. In addition, in ther- 
mal equilibrium the i?, A and K components of the 
bosonic and fermionic propagators are related by the 
fluctuation-dissipation theorem (FDT). Below we demon- 
strate that our approximation, Eqs. ( pTj) and (^8|) , for the 



Green function satisfies these requirements. For this pur- 
pose it is convenient to rewrite identically Eq. (^^ in the 
following form 



7>«) 



>«)g>«) 



where 



B'^{t)^B^{t)=B>{t)-B<{t), 
B^(t)=B>{t) + B<{t). 



(71) 



(72a) 
(72b) 



The > and < components of the fermionic Green func- 
tions are related to the R, A and K in the same manner. 
^From Eqs. ( [72|) and Eqs. (^) we obtain 



B>i<) (CO) = Idt e^- exp U / e^"'* 5:(V«(g, co') - V^{q, u')) ( coth ^ ± 1 



2tt 



2r 



(73) 



According to the FDT the equilibrium bosonic and 
fermionic Green functions in the frequency representa- 
tion satisfy the following relations 



B>{uj) =exp{tj/r}6<(tj), 
g>{e)=-exp{e/T}g<{e) 



(74a) 
(74b) 



It is not difficult to see that if any pair of bosonic Green 
functions [t) and B'^{t) satisfies Eq. (74a) then for 
any analytic function /(z) the pair f^{t) = f{B^{t)) 
and /^(t) = f{B^{t)) also satisfies it. Indeed, 



/>«)(c.) 



( / 



)>(<) /, r\-iu)'t 



(75) 



Expanding / on the r.h.s. of this equation in the Tay- 
lor series and performing the t integration, we see that 
in each order of the expansion f^{uj) = exp{u>/T)f^{uj). 
One can also check that if C/^^^^ and S^^^-* satisfy the 
FDT, Eq. (H), then so do the functions ^>«) defined 
as 



^>«)(t) =g>«)(i)S>«)(i), 



(76) 



9 



Noting that the arguments in the exponential in Eq. j?^ ) 
obviously satisfy the FDT, Eq. (74a), we conclude that 
B-^^'^\ll>), Eq. (|73|), and the approximate Green function, 
Eq. (|67|), satisfies it as well. 

To establish the positive definiteness of the tunneling 
density of states, Eq. (|70|), we first show that B''^^\u) in 
Eq. (|7^) is positive definite. Indeed, exp[— jV^(O)] is real 
as can be seen from Eqs. (|65|). It is also not difficult to see 
that each Fourier component of the argument of the ex- 
ponential in ( [73| ) is positive definite. All the coefficients 
in the Taylor series of the exponential are positive and, 
since the Fourier transform of a product is given by the 
convolution of positively defined Fourier transforms, we 
conclude that the l.h.s. of (|73) is positive definite. We 



next use Eqs. (72a) and (74a) to rewrite the tunneling 
density of states as 



(77) 



where 5^(e) is given by Eq. (|7l|). Since Af > we im- 
mediately see from Eq. (|7l|) that the tunneling density of 
states is positive. 

In equilibrium, it is convenient to write the DOS 
through the Keldysh Green function employing the FDT 



tanhe/(2T) 



(78) 



As was proven above, Eqs. ( [77| ) and (jTl) are equivalent. 
One can then express B^{t) through B^^\t), where the 
latter are conveniently rewritten as 



B>«)(t) = -exp 




(coth^(l 



2T' 



COS, Lut) ± i smut] 3 V^'{q 



(79) 



Expanding this expression to the first order in the in- 
teraction, V, and substituting into Eq. ( |78| ) one recovers 
the Altshuler-Aronov result for the zero-bias anomaly 
p6[ . This perturbative result corresponds to the diagram 
drawn in Fig.^. 




FIG. y Lowest order interaction correction to the single par- 
ticle Green function. The wavy line here denotes the RPA- 
screened Coulomb interaction. The impurity-dressed single 
particle Green functions are depicted by solid lines and the 
double dashed lines represent diffusons. 



We shall restrict ourselves to the analysis of the non- 
perturbative result, Eq. (^, (|7|), only at T = 0. Noting 
that for T — 0, Ft — {int)'^, one obtains 




In the two dimensional case Eq. ( |65a| ) with Vo{q) — 
27re^/(7 leads to 

oo 

/^3^V'(,.„)('-— ') (81, 
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_ 1 ( Int/rlntTUJ^ + 2^\ntLjQ 
SiT^g \ irlntLUo ' 

where g = vD is the conductance; = Dk^ and 
K = 27re^i^ is the inverse screening radius; 7 — .577 ... is 
the Euler constant. Since we have neglected the fluctua- 
tions of Q, we have missed corrections of order 3"^ Ini/r 
(in d = 2), therefore Eq. ( |80| ) can only be trusted for e 
not too small, such that (STr^g)"^ ln(er)~"'^ ^ 1. How- 
ever In^ t/r terms have been accounted for correctly by 
our procedure. If in addition g~'^\n.ujQT <C 1, the time 
integral in Eq. ( pO| ) may be performed by the stationary 
point method, resulting in 

Ke) = ^exp|-^ln(|e|r)-iln™oV|e|| . (82) 

Theoretically, however, g^^lnuj^r, need not be small. 
In that case the stationary point integration should be 
somewhat modified and terms oc Intwg should be re- 
tained. 

We have achieved a non-perturbative resummation 
of anomalously divergent, cx In^ er, terms for a single 
particle Green function. The non-perturbative expres- 
sion for the DOS, essentially arises from the gauge non- 
invariance of the single particle Green function. The 
calculations above are in essence the "Debaye-Waller" 
factor 1^ due to almost pure gauge fluctuations of elec- 
tric fields, cf. Eq. (|62|). Gauge invariant characteristics 
(such as e.g. conductivity) do not carry phase factors 
and therefore are not affected by the interactions on this 
level of accuracy (fluctuations of Q should be retained) . 
In perturbation theory this fact is reflected by the can- 
cellation of diagrams without diffusons (apart from those 
which renormalize vertices) [p6[. In our formulation such 
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terms never appear since the phase factors cancel along 
any closed loop diagram. 

The gauge physics of anomalous corrections to the 
DOS was first realized by Finkel'stein [ p"2|j3l| ], who ob- 
tained a non-perturbative result similar to ours. Nazarov 
and later Levitov and Shytov obtained the same 
result (in imaginary time) by semiclassical reasoning. 
Kopietz recently reinstated it, stressing the role of 
phase fluctuations. The analogous expression for the 
zero-dimensional case has also been known for some time 
p^-p6[. We believe that we provide its first consistent 
derivation using the cr-model. Unlike the previous ap- 
proaches, the Keldysh technique provides the answer di- 
rectly in real time and finite temperature. This enables 
us to circumvent the tedious analytical continuation pro- 
cedure. 



B. Shot noise 



In this subsection we shall use the results obtained in 



subsection [IV A| to calculate the power spectrum of cur- 
rent noise through a tunneling contact between a clean 
metal and a dirty metal film. The power spectrum of 
current noise is given by 



Siij) = / dte^"*((/(<)/(0) +/(0)/(t)» 



(83) 



Here the current operator in the tunneling approximation 
is given by 1(1) = iTa\t, ro)bit, tq) - iT*b^ [t, ro)a(t, ro), 
and 6(t, ro) and a{t,ro) are electron annihilation opera- 
tors at the position of the contact in the dirty film and 
in the clean metal respectively. Below all the fermion 
operators and the Green functions are taken at the point 
of the tunneling contact tq, and we omit the position ar- 
gument for brevity. Using the expression for the current 
operator we can rewrite Eq. (k3) as 



s{l.) = |Tp/dte-* g<(-t)g>it) + g>{t)g<{^t) + g<it)g>{-t) + g>i-t)g<it) 



(84) 



where ga and C/h are Green functions for the clean metal and for the dirty film respectively. We assume that the voltage 
V is applied across the contact. To the lowest order in the tunneling matrix element, the Green functions under these 
conditions are equilibrium, except that the chemical potentials in the two metals differ by eV. Therefore in the lowest 
order in the tunneling amplitude we can express the power spectrum of current noise through the equilibrium Green 
functions. Expressing them through DOS with the aid of the FDT and utilizing the fact that for the clean metal 
DOS, i^a, is independent of energy e one obtains 



S{uj) = 2Triya\T\^ devbie) n{e)[2 ~ n{e + uj - eV) - nie - lu - eV)] + [1 - n(e)] (n(e -|- w - eU) + n(e - 



eV)) 



(85) 



r 



where n(e) = [1 -I- exp(e/T)]~^ is the Fermi function. 

Setting y = in Eq. ( |85| ) we obtain the power 
spectrum of the equilibrium current noise in the con- 
tact So{uj). The excess noise is given by the difference 
dS{uj) = S{u;) — So{lli). The noise power is a symmetric 
function of frequency, and at zero temperature reduces 
to 



6S{uj > 0) = 2Triya\T\ 



Lj+eV 



-eV\ 



(86) 

At zero frequency the shot noise is proportional to the 
total current. This is natural, since in the lowest order in 
the tunneling amplitude the electrons pass through the 
contact extremely rarely and, therefore, can be viewed as 
noninteracting. The role of interactions reduces to mod- 
ification of the density of the available states. The cusp 
present at zero temperature in the noise power spectrum 



for noninteracting electrons at a; = and lo — eV is 
washed out because DOS vanishes at e = 0. 



C. External fields and auxiliary sources 

In some sense our previous manipulations leading to 
Eqs. dsi]), ( p2| ) were no more than a complicated rep- 
resentation of unity. This is so since, according to the 
basic idea of the Keldysh technique, the partition func- 
tion, Z, is identically equal to unity. To make the entire 
construction meaningful one should introduce auxiliary 
source fields, which enable one to compute various ob- 
servables. We shall do it in parallel with introducing ex- 
ternal classical fields. Since we shall mostly discuss the 
conductivity, we'll use the vector potential A(r, i) as an 
example Other fields (e.g. the scalar potential) may 
be introduced in a similar way. We introduce a doublet 
in the rotated Keldysh frame 
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A(r, t) 



ai(r, t) 
a2(r, t) 



(87) 



which is related by the usual transformation (cf. 
Eq. (p^)) with the two fields aa{r,t) residing on the two 
branches of the time contour. The vector potentials en- 
ters the fermionic Hamiltonian through the long spatial 
derivatives, Vr ^ Vr + ikaj". The classical external 
vector potential is the same on the two branches of the 
contour and hence it is described by the symmetric com- 
ponent, ai(r, t), only, whereas a2 = 0. In this case the 
generating function still equals to unity 

Z[ai,a2 -0] = 1. (88) 

To obtain a non-trivial generating function one has to 
introduce a quantum component of the source field, 
a2 (r, t) . This component does not have a classical mean- 
ing and thus has to be nullified at the end of the calcu- 
lations. Its presence however is essential for generating 
observables. One can easily check that the current den- 
sity defined as 



1=1,2 




Here we have omitted the vector indices of a using its lon- 
gitudinal character. In general, any response function is 
given by the (2,1) component of the appropriate bosonic 
correlator. The structure of the theory guarantees that 
this is a retarded function (cf. e.g. Eq. (p4[)). 

In the presence of an external vector potential, A, the 
trial saddle point, Eq. (|30|), is shifted. Noting that A en- 
ters the action always in the combination V K + A, one 
finds that the condition for the optimal K is given by 
Eq. d^), with the substitution Dq'^K Dq^K + iD{qA) 
p9| . Solving this equation in a manner it was done in sec- 
tion [II A , one obtains that Eq. ( ^l| ) should be modified 



V-\q,Lu)K{q,u) = njl$(g,t^) - iDai{qA{q,,j)) , (93) 

where bosonic propagators uj) and = —iuj'D{q = 
Q,uj) are defined by Eqs. (^) and (^4|). In solving 
Eq. (37) with the external field we still assumed that 
the fermionic distribution function is the equilibrium one. 
This is a legitimate procedure in linear response. The 
generalization to the non-equilibrium case is discussed 
in section VI . After disregarding the massive modes and 
expanding the logarithm, one obtains Eq. (^8|), with the 
long derivative modified as 



is given by [R8| 



./ N e 5Z\A] 
2i dai2(r,t) 



(90) 



32=0 



We restrict ourselves to the case of longitudinal vector 
potentials only. Taking into account the fact that the 
external electric field is given by — iu;ai(q, aj)/e, one ob- 
tains that the linear response conductivity is given by the 
retarded component of the current-current correlator 



where 



a{q,uj) 



(9,w) = — 



2i (Sa^(q, a;)(5aa(— q, — oj) 



(91) 



(92) 



A=0 



drQ = VrQ + i[(Vfc„ + a„)7", Q] 



(94) 



and K given by Eq. (|93|). Since 7^ = 1, any static exter- 
nal field, ai(r), appears to be decoupled from Q. This 
refiects the fact that diffusons are not coupled to a static 
magnetic field. On the other hand, even space and time 
independent quantum component, a2, is coupled to Q. A 
little algebra shows that 

VK + A= -iqVn-^ ($ + {qA)Lu/q^) ; (95a) 
$ + iujK= Dq^Vai ($ + {qA)uj/q^) . (95b) 



With these expressions and the long derivative given by 
Eq. (|9^), one can rearrange Eq. ( |48| ) to obtain the average 
generating function in the following form 



{Z[A]) = y"l?$cxp{2Tr{$^V'^i* + [^+M)^/9']^^o[* + (9^)^/9']}} j VQei^^\^^iSi[Q,V K + A 



(96) 



r 



where the action Si, I — 0,1,2 is given by Eqs. (|5^ ) 
and the bare polarization operator, Po{q,uj), is given 
by Eq. (^5|). By virtue of Eqs. ( p5|) the entire action 
is expressible through the combination $ + {qA)uj/q^, 
which is proportional to the (gauge invariant) electric 



field V<I> 4- dtA. This fact immediately guarantees that 
the continuity equation is satisfied to all orders in the 
perturbation theory. Indeed, one could introduce the ex- 
ternal scalar potential, (p, which enters the action always 
as $ -|- (/9 (apart from the bare interaction term, Vq)- 
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Then, due to the fact that Z = Z[ip + {qA)uj/q% the 
charge density, p = {2i)~^5Z / 5ip2-, and current density, 
j = —{2i)~^5Z/5sL2, has to be related by 



(Vj)+9tp = 0. 



(97) 



The corresponding variation with respect to the classical 
components, (pi and ai, guarantees continuity at each 
branch of the contour separately even in the presence of 
non-zero auxiliary quantum fields. As a result of con- 
tinuity the exact relations between current-current and 
density-density correlators holds 



(98) 



At the saddle point, Q = A, one has Si[A, VK+A] = 0. 
Thus, neglecting the fluctuations of Q, one obtains for the 
RPA generating function 



exp 



{Zrpa[A]) 

T>^ exp |iTr{$^y"i$ + 2$^^Po^} 



(99) 



Performing finally the Gaussian integration, one finds 

{Zrpa[A]) = exp |zTr{^^^P«p^A}| , (100) 

where PppAiq,^^) is the RPA screened density-density 
(polarization) correlator, which is given by 



PRPAiq.io) = Po - Po[Vo-' + Por'Po - [Po 



and has the structure of a bosonic correlator 



Vo] 



Piq,to) 



P^(9,c^) 
P«(g,a;) P^(g,c^) 



(101) 



(102) 



The up' jq^ factor in Eq. (|lO0| ) reflects the relation be- 
tween the density-density and the longitudinal compo- 
nent of the current-current correlators, Eq. (|9^). The 
fact that the (1,1) component of P (as well as of any other 
bosonic correlator) vanishes is a manifestation of the nor- 
malization condition, Eq. (^). Employing Eq. ( pl[ ) one 
obtains for the conductivity in the RPA 



(7RPAiq,uj) = e^vD 



Dq''{l^vVQ{q))-iuj 



(103) 



One is usually interested in the irreducible part of the 
density-density (or current-current) correlators, which 
describes the linear response to the total or internal field 
and not to the external field as discussed above. The re- 
lation between the irreducible part, P^^r, and the total P 
is exactly the same as between the bare, Pq, and Prpa, 



Eq. (101) 



(104) 



Shifting the integration variable $ ^ $ — {qA)u!/q'^ in 
Eq. (^) and differentiating twice with respect to A one 
obtains an exact relation for the polarization operator 

Piq, Lu) = V' + 2*K,"'(<i'(<7, co)'S>^i-q, -cc)) , 

(105) 

where ($$^) is an exact propagator (averaged wit h rc - 
spect to the full action, Eq. (96)). Employing Eq. (104) 
one finds 



P.rr-'^{{^'f'^)y'-V,-' 



(106) 



If one is interested in the response to a uniform exter- 
nal field, q = 0, the expressions may be further simplified. 
Noting that for the Coulomb interaction Vf^^{q = 0) = 
and employing Eq. ( psf ) and relation between <i> and K, 
Eq. (pl|), one obtains 



Y.,„{q = 0,Lu) 



(107) 



Unlike in the calculations of the single-particle Green 
function, only \7 K and never K itself appears in cal- 
culations of gauge invariant quantities. This allows 
one to consider a universal limit of strong interactions 
^o'^iq) ~^ t^i^ limit it is convenient to change the 

integration variable from $ to V K (although formally it 
is a vector it has only a longitudinal component and hence 
number of variables is conserved). In the new variables 
the Gaussian weight is given by iTiiS/ q~'^V^^'^ K} , 
where V is defined by Eqs. (64)-(|65|). In the universal 
limit one has 



V-\q,uj)~^-iyDq^V-\q,u;), 



(108) 



where the diffusion propagator V is defined by Eqs. (|4^), 
(^). Finally, one obtains for the action in terms of VK 

2 

I- I- 'yiSilQ^VK] 

{Z) = VWK e— cTr{VK-p-iVK} /^g ^to 

(109) 



Eqs. (107) and (10£) constitute a complete framework 
for calculations of gauge-invariant response functions. 
Neglecting Q-fluctuations one finds Eirr(9 = 0,w) = 
vDiujIiy , which leads to the Drude conductivity, a =^ 
e^vD. Fluctuations of Q and VK lead to weak- 
localization and interaction corrections. Note that un- 



like in the case of the DOS (section IV A), fluctuations 
of VK alone, with Q = A, do not lead to any corrections 
to linear response. This is a direct consequence of gauge 
invariance of linear response functions. Only combined 
fluctuations of V-fiT and Q, discussed in the next section, 
renormalize the Drude conductivity. 
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V. FLUCTUATION EFFECTS 



A. Q— matrix parameterization 



As was discussed in section [II B the massless fluctua- 



tions of the Q-matrix can be parameterized as 



Q = exp{-VF/2}Aexp{M^/2} , 



(110) 



where 



saddle point. As it will be apparent later, the physical 
meaning of w is a deviation of the fermionic distribution 
function, F, from its stationary value. At the same time, 
w has no classical interpretation. To a large extent it 
plays the role of the quantum counterpart of w, which 
appears only as the internal line in the diagrams. 

One may expand now the action, Eqs. (p2|), in powers 
of w and w. The expansion of the non-interacting action, 
jS'o[(5] starts from the second order, which has a familiar 
diffusive structure 



WA + AW = 0. 



(Ill) 



-We 



We 



(113) 



Employing Eq. ( p4| ) , one obtains that the general form of 
W, which satisfies the condition (111) is 



W = 



1 F 

-1 

Fw - 
—w 



w 
w 

-w + FwF 
-wF 



1 F 

-1 



(112) 



As a result the bare propagator of the Q-matrix fluctu- 
ations is given by 



{We^eAq)Wc3cA-<l))w = " ITT ^'^'^^'^"^ 



TTiy Dq^ + i(ei - 62) 

2^eie3^e2e4 r^At \ 
D {q,6i-62). 



(114) 



where w and w are arbitrary Hermitian matrices in the 
time space. Below we shall thus understand the func- 
tional integration over Q as integration over Hermitian 
w and w. Notice that Q itself (as well as the Green func- 
tion G) appears to be non-Hermitian. It means that the 
"contour" of integration in the Q space is deformed from 
being pure Hermitian to pass through the non-Hermitian 



The higher order terms describe non-linear interactions 
of diffusons with the vertices having the structure of 
Hikami boxes. One can easily work out this expansion 
in the Keldysh formalism. We shall not do it here, since 
our main focus is on the interaction effects. Substitut- 
ing = [A, W]/2 into iSi[Q., VK], one obtains in the 
first order in W 



iS^^\w,VK] 



-^Tr{[i^V^A:„(A7"A-7") + 



iwfc„)(7"A- A7")] W] 



In equilibrium iS\^'^\W,\I K] = 0. Indeed the r.h.s. of Eq. (|115|) coincides with equation 



(115) 



which was used to 



determine the K functional. In equilibrium we were able to solve Eq. ( pT] ) by an appropriate choice of K . This was 
precisely the motivation behind looking for the saddle point for each realization of the Hubbard-Stratonovich fleld: 
to cancel terms linear in W . Since we could not find the exact saddle point, such terms do appear, however, only in 
the second order in VK. For ^S2[Q^'^ K] part of the action one obtains 



iS'i\w, VK] = ^VA;„(ei ~ e2)Tr { [7"A,,7'5A,3 - A,,7"A,,7''] W} Vkp{62 - £3) 



■KvD 



nvDTr {Vif^(ei - 62) [M, 



M-e,.,^e,e 



J VK{62 - ea)} , 



where we have introduced two vertex matrices in the bosonic Keldysh space 
1 \ / 2Fe.- F„ - F, 



M. 



^+Fe,Fe, 



1 -2Fe, J ' ^"^1^^^^ ~ \-l~Fe,Fe,+ 2fJf,, F^, + 



IF F 
IF F F 



(116) 



(117) 



The fact that (M, 



w\l,l 



is a manifestation of the normalization condition, Z = 1. Indeed, this matrix element 



connects only the classical components of W and K fields, which alone can not change the normalization. Being 
averaged over fluctuations of Vif with the action Eq. (109), iS'2^''[VF, ViT] gives 



Si'\W,VK])vK = 27nWee [iFe+^ - Fe)V'^ij) - (1 - Fe+^ Fe){V'^' (lj) - (io))] 



(118) 



There is no term proportional to the classic component, w. In equilibrium the r.h.s. of Eq. (118) is obviously zero. 
Out of equilibrium, it is this term which is responsible for the standard collision integral, see section VI. As we shall 
see in the next section, fluctuations described by [W, VK] are responsible for the Altshuler-Aronov corrections 
to conductivity. For completeness we write also the second order expansion of iSi[W,\i'K] 
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(2) 

iSl '[W,VK] = ii^D V/ci(ei - e2)(Vwe2£3'^«3«i - ^f'€2€3^'^«3ei ) + 



(119) 



B. Altshuler— Aronov corrections 



Restricting oneself to the lowest non-vanishing terms in the expansion over W, Eqs. (113) and (116), one obtains 
a Gaussian theory with respect to the W fluctuations. After integrating out these fluctuations employing Eq. (114), 
one ends up with the action for the WK field only 



(120) 



27T,.D'Tr {V<_,,(r)M™,^,3Vif,,^,3(r)} D^{r - r' , eg - ei)Tr {V<_,,(r')Af-VX,,_,, (r')} , 



This way the (Vif )■* effective vertex is generated. Perturbatively the {WK)* interaction term may be treated by pairing 
two fields, say V {V KV K^)V K . This results in a renormalization of the bare correlator, 'D~^ . Only pairing of 
WK fields in different spatial points leads to non-vanishing corrections, see Fig. |[ There are four different ways one 
can pair (VfcQ,(r)V/c^(r')). Taking into account all these four possibilities and integrating over an intermediate energy 
one obtains correction for e.g. retarded component of the {{V KV K'^))"^ correlator 



(121) 



q' ,uj' 



r 



where B^^ is defined by Eq. ( pq ) ; d is dimensionality. Ob- 
viously, the correction preserves the retarded character of 
the corresponding component. In equilibrium, the cor- 
rection to the Keldysh component obeys the fluctuation- 
dissipation relation 



[sv-Y = coth^[[sv-r-[sv-Y] 



(122) 




FIG. ^Lowest order self-energy diagram for {S/ KS/ K"^} prop- 
agator. The zigzag lines represent the bare {\/K\/K'^) prop- 
agators, the parallel solid lines denote the (WW) propagator 
and the open circles with two zigzag and two straight lines 
emanating from them represents the VK^WVK vertices. 



Employing Eqs. (|9l|), (107) one obtains for the correc- 
tion to the q — conductivity 



Ae^D 



[{lu' - uj)B^' 



D'^W, ^ + ^')[D{q', u;'Y X (123) 



In the low frequency limit this reduces to the familiar 
expression [Eg 



oo 

5<T = i — - / dijj —— (oj coth — I 



(124) 



where <Jd = e^vD. 





+ 



FIG. ^ Lowest order diagrams for the interaction correction 
to conductivity. Their sum is equivalent to the diagram in 
Fig. ^ in the present formalism. 

Note that this expression is given by the sum of dia- 
grams drawn in Fig. ^. The other diagrams which are 
presented in Fig. |^ add up to zero. They represent the 
purely phase correction to the single particle Green func- 
tion and therefore do not enter the expression for the 
conductivity. In the present formalism these diagrams 
do not appear at all. 
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In two dimensions expression (124) leads to the loga- 
rithmically divergent negative correction to the conduc- 
tivity (or conductance) 



+ 





+ 



FIG. IE] Diagrams for the interaction corrections to conduc- 
tivity which add up to zero. These diagrams never appear in 
our formalism. 



6a_ 
a 



9 



■InTr^ 



(125) 



where the elastic mean free time, r*^', enters as an upper 
cutoff in the integral over frequency. 

To handle this divergence one may try to setup a 
self-consistent mean-field treatment of the (WK)^ non- 
linearity. To this end let us put {V KV K^) propaga- 
tor on Fig. ^ to be a dr essed one, 2?, where f)^^ — 
V-^ + SV"^. Then Eq. (|l2ll) may be rewritten as a 
closed non-linear equation for e.g. retarded component 
of the propagator, [2?(g.w)]^, 



dv 



J2 D'^iq + q',0J + u')[V{q', [{u:' - u;)B^,^^ - J B^,] 



q' ,cj' 



[V{q,L0) 



(126) 



The frequency dependent conductivity is then given by 

([p(g=o,..)n-i 



cr(w) = e vD ■ 



(127) 



One may easily check that in the one-loop approximation 
there are no other corrections to conductivity. Indeed, 

the possible corrections like {{iS^'^\W,V K]^ )^ = 

and {182 \W,'\/ K])w = 0. They vanish since they in- 
clude energy integration of purely retarded or advanced 
functions. Being expanded to higher orders in W these 
terms yield weak-localization corrections. 



VI. KINETIC EQUATION 

The aim of this section is to demonstrate how the ki- 
netic equation for the distribution function F appears 
naturally in the framework of the Keldysh formulation. 
The kinetic equation is nothing but the saddle point 
equation for the effective action on the Q-matrix pO| . 
In the case of interacting electrons it is obtained by inte- 
grating out the K (or equivalently <I>) degrees of freedom. 
Consider the partition function, Eq. (21a), with the ac- 
tion, S{Q,<^), given by Eq. (H). Let us perform the <I> 
integration first. As a result we obtain for the average 
partition function 



{Z) = J VQe'^'-ff^'^^ , (128a) 
iSefflQ] = In J p$e'Tr{$-i/-i.i*}+.s[<5,*] (128b) 



Since the the action, S{Q, <&), Eq. (^), is quadratic in <& 
(given the linear relation between K and $) the integra- 
tion in the last expression can be carried out explicitly. 
We find it more convenient, however, to proceed with the 



expression (128b). To obtain a non-trivial kinetic theory 
one may assume the presence of classical external fields, 
like e.g. scalar or vector potentials. These fields can be 
introduced in the action Eq. (^) in the way it was done 
in section 



IV C 



We shall look now for the saddle point equation for Q 



SS, 



effl 



SQ 



(129) 



obtained under the condition — 1. Let us reiterate 
again the logic of the entire procedure. After averaging 
over disorder and introducing the Q-matrix, we found 
that the low-energy degrees of freedom are described by 
the Q-matrices given by Eq. (|3l| ) with — 1. We then 
restrict ourselves to this massless manifold and look for 
a realization of Q which extremize the effective action. 
The latter is obtained by integrating out the photon fields 
originating from c-e interactions. Without any external 
fields (and/or non-trivial boundary conditions) such an 
extremal Q is simply given by A, Eq. (|2^), with the equi- 
librium F function, Eq. (26). If external fields (and/or 
non-trivial boundary conditions) are present, the sta- 
tionary Q may deviate from A, still being on the massless 
manifold, = 1. The stationary point is to be found 
by solving Eq. (129), which turns out to be precisely the 



kinetic equation with the collision integral term. 

Before proceeding along these lines, let us comment 
on the relation between the phase K, introduced in sec- 
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tion III A and the Hubbard-Stratonovich field, The 



procedure of section 111 A was based on the property of 
the equilibrium distribution described by Eq. (|40|). We 
need to generalize it for non-equihbrium situations. To 
this end we note that the equation for the quantum 
component, k2{r,t), Eq. (p8|), does not contain distri- 
bution function and remains valid for a non-equilibrium 
case. The equation for the classical component ki{r,t), 
Eq. (|9|), can not be satisfied identically out of equihb- 
rium. Thus the choice of fci (r, t) allows for a certain 
arbitrariness. However, as we shall see below, this arbi- 
trariness does not affect the form of the kinetic (saddle 
point) equation. It would manifest itself in a calculation 
of fluctuation corrections (cf. section ) to the non- 
equilibrium saddle point result. We shall not attempt 
this task here. For our purposes it is sufficient to keep 
the definition of K{r, t) given by Eq. (|l|) (or Eq. (||) if 
external fields are present). The equilibrium bosonic dis- 
tribution function, used in the definition of the Keldysh 
component of the propagator 2?(g,w), Eq. (^3|b), does 
not show up in the kinetic equation. 

Employing Eq. ( 128b| ), we rewrite the saddle point 



equation (129) as 



. SQ 



where 



(•••). 



(130) 



(131) 



Here Q is a self-consistent saddle point solution of 
Eq. ( |130| ). Performing variation of the action 5'((5, $) 
given by Eq. (|8|) under the condition = 1, one ob- 
tains 



(132) 



where Q — 1. This equation is analogous to the kinetic 
equation in the semiclassical theory of disordered super- 
conductors |4^,^^. We have derived it here for the case 
of a normal interacting metal. 



We shall seek the solution of Eq. (132) in the classical 

form, e.g. obeying the condition Q^_^ = 0. A non-zero 

quantum component at the saddle point would violate 

- 2 

causality. Provided Q^^ = and Q = 1 are satisfied the 
saddle point solution assumes the form 



af^(e-e') 2F,,Ar) 
-lf(5(e-e') 



(133) 



where F^^^i^r) is a non-stationary distribution function. 
Assuming that the saddle point has the form given by 
Eq. ( |133| ), one can easily check that the exponent in the 
<i> averaging, Eq. (131), does not contain linear terms in <I> 
(or VK). Indeed, the terms proportional to V/ci vanish 
identically, which is a manifestation of the normalization 
condition, Z = 1. From another hand, terms propor- 
tional to Vfc2 are reduced to the full gradient (the fact 
that there is no ambiguity in the choice of ^2 is important 
here) and thus also vanish upon the spatial integration. 
As a result the t erm s linear in $ (or WK) in the saddle 
point equation (|l32| ) do not survive the $ integration. 
Therefore Eq. (|l32| ) may be reduced to 



DVriQVrQ) 



■ i[e,Q] = D{[VkarQyki3j^,Q])<s> . 

(134) 



The r.h.s. of this equation contains the collision integral 
term along with the coUisionless renormalization of the 
kinetic part. To evaluate it one needs to know the propa- 
gator {Wka{r, t)Wk/s{r' , i'))$ at r = r', averaged over the 
non-equilibrium action, Eq. ( 131 ). To follow the same 
notations as for the equilibrium case we shall denote this 
propagator as 



(Vfc„(r,t)Vfc^(r',t'))* 



'^K^'iry). (135) 



The form of the saddle point, Q given by Eq. (133) guar- 
antees that Vt^t' has the standard retarded/advanced 
structure of a Keldysh propagator. Employing Eq. (|2^), 
one finds that the only non-zero matrix component of 
Eq. (134) is its Keldysh (1,2) component. The corre- 
sponding equation for the distribution function Ft^t' (^) 
takes the following form 



DVlFt 



t.t' - [dt + dt')Ft.t, = - 

V 



Ft_ 



V 



K 



(136) 



Here all F functions and propagators V are to be taken at the same spatial point; integration over ti is assumed in 
the last term on the r.h.s. Note that the l.h.s. of this equation is a linear diffusion operator acting on Ft^t'ir)- The 
subsequent calculations are significantly simplified by passing to the Wigner representation. 



F,(r,r) ^ j j dtdt'Ft,t'{r)e'''-'-*'^S (^i 



t + t' 



(137) 



Furthermore we shall assume that i^e(r, t) is a slow function of r on the scale 1/T (or any other inverse characteristic 
scale of energy, e). With this assumption Eq. (|136D may be rewritten as 
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DVIF,{t) - drF,{T) - 9,7^,(r)5eF,(T) + d,'R,{T)drF,{T) 



(138) 



where 



J^.-c.(r-,T) 



(139) 



The r.h.s. of Eq. ( |l3q ) represents the coUision integral, of. Eq. (118). If equihbriuni relation, Eq. ([I3|b), between 
Keldysh and retarded and advanced components of T) holds, then the equilibrium distribution function, Eq. ([2^), 
nullifies the collision integral. Therefore Eq. (138) is satisfied in the thermal equilibrium. This, in fact, provides 
justification for our previous use of A with the equilibrium i^-function as the saddle point. Indeed, without interactions 
(and hence without collision integral) any stationary function F^^ satisfies the saddle point equation. It is the relaxation 
processes due to e-e interactions that render the equilibrium solution unique. The terms which contain real part of 
the self-energy, 7?.e(r, r), lead to a coUisionless renormalization of the kinetic part, see section [VI B . 



To proceed further we need an explicit form of the non-equilibrium propagator, I?i^(r, r, r). We shall evaluate it 
in the universal limit of strong interactions, V^^ 0. Substituting the saddle point Q given by Eq. (133) into the 
action S[Q,^], Eq. (Eq), and performing the Gaussian integration one finds (cf. Eq. (pfl) ) 



^DVlaf + dry f E {7"a+^ (r, rh^Q^_^ (r, r) - 7-7^} 



(140) 



The term with originates from the term + pcji in the action Eq. (|4^). (It is easy to check that the ambiguity 
in the choice of ki, mentioned above, disappears upon the calculation of this term by the symmetry reason.) The local 
in space term in Eq. ( |140| ) originates from DTT{drQ)'^ ■ Assuming that any distortion of the equilibrium distribution 
is limited to a vicinity of the Fermi energy, e.g. Ff:^±oo{r,T) —> ±1, one finds 







-L'V^ + iuSr 



— -DV^ — itoSry —2iujSr,r'Buj{r, r) 
By definition, the non-equilibrium bosonic distribution function is given by (cf. Eq. (|5^ ) ) 



(141) 



OG 

BUr, r) = ^ Jde[l- F,+ ^{r, r)F,_^ (r, r) 



(142) 



According to Eq. (141) the retarded and advanced components of V are not modified with respect to their equilibrium 
value, Eq. (^3|a). As a result, Pf^^"^' (r, r', t) = Pf^'"*' (r— r') even in a non-equilibrium situation. Inverting the operator 
on the r.h.s. of Eq. (141), one finds for the Keldysh component at coinciding spatial points 



VS{r-r')BUr',r)V^ir' -r) + 



(143) 



2i 



VSir - r')drBjr',T)d^V^{r' - r) - d^V^ir ~ r')drB^{r\ T)V^{r' - r) 



^From now on we shall retain only the first term in this expression, which is dominant due to the assumed slowness 
of the temporal variations of F^ (r) . If in addition (r, r) changes slowly on the spatial scale L^^ = ^JdJuj, where 
LD T, the expression for the Keldysh component acquires the quasi-equilibrium form 



r, r) ^ i?.(r, T)Y,[V^{q) - V^{q)] 



(144) 



One can calculate gradient corrections to this expression, which lead to a non-local collision integral. Usually such 
corrections may be safely neglected. Finally in this hydrodynamic regime the kinetic equation takes the form 
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DVlF.ir) - [1 - dMT)]drF,iT) - dMT)d,F,{T) 



(145) 



E 



BUrm-Ur) - F^ir)) + (1 - F,_Ur)F,{r)) 



with 



TZ,{r,r) = -J2[^VS{q)]F,.Ur,r) 



(146) 



and Vl^iq) = (Dg^ _ iuj)^^. 



A. Collision integral and relaxation time 



Using the conventional fermion distribution function, n^{r, r) — {1 — Ff^{r, t))/2, one can rewrite the colhsion integral 
in the usual form with "out" and "in" relaxation terms. Indeed, employing Eqs. (|57|), (|5^), one identically rewrites 
the r.h.s. of Eq. (|145[) as 



dLuae q 



,{l - ne/)(l - rie-t^) - ne/n£_^^(l - n^){l - n^'-^) 



(147) 



This is precisely the collision term derived by Altshuler and Altshuler and Aronov two decades ago. One can 
linearize this expression around the equilibrium distribution by the substitution 

F,{r,T)^F!'i-w,{r,T)/2 (148) 

and keeping linear terms in Wt{r, r). This way one derives the familiar results for the relaxation rates Q,^,^. We 
shall not repeat this procedure here. Instead we shall demonstrate how this quantities may be extracted directly from 
the effective action. To this end we note that the kinetic equation ( |13^ ) maybe written as {2Tny)^^SiSeff /Swe{r, r) = 0. 
(As usual, an observable is generated by differentiation with respect to a quantum component.) Thus the linearized 
version of the kinetic equation is just —{nv)~^6'^iSef f /SwSw\Q^f^. According to Eqs. ( 128b| ) and (131) 



SwSw 



SwSw 



SiS[Q, $] SiS[Q, $] 



SiS[Q, $] \ / SiS[Q, $] 



Sw 



Sw 



Sw 



Sw 



(149) 



where all the variational derivatives are calculated at Q = Q = A . The last term in this expression vanishes 
identically, since A is obviously a solution of the kinetic equation (130). The first term originates from the expansion 
of {iS2[W,W K])yK , Eq- ( |52c|) , to the leading order in w and w. After a little algebra one obtains (there are no terms 
with WW in equilibrium) 



) P«(B:^ + J - V^iBZ'^ + Fl 



(150) 



where e± = e ± Q./2. Equivalently this expression can be obtained by variation of the F^(r, t) functions in Eq. ( |136[ ). 
The terms with We^,c_ and We^^-u},t_-u} represent "out" and "in" relaxation processes respectively. Their condensed 
diagrammatic representation is given in Fig. |^a,b. The full set of corresponding original diagrams may be found e.g. 
in Ref. E^. Restricting ourselves to the diagonal fluctuations, = 0, only, we obtain for the "out" relaxation rate 



1 



E 

ujq 



V 



At T = in two dimensions this leads to the familiar result [fe6| 



,0; e ~ uj 

coth h tanh 

2T 2T 



(151) 



1 

^d=2 



(e) 47r(7 ' 



(152) 
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where g = vD. Expanding the r.h.s. of Eq. ( |150| ) in a small fJ, one can also recover the collisionless terms in the l.h.s. 
of Eq. dli^ ). 

Finally we concentrate on the second term in the r.h.s. of Eq. (149). This term corresponds to the variation of 



T)^^, in Eq. (136) over a deviation from its equilibrium value (or equivalently variation of B^{r,T) in Eq. (145) ). 
Its condensed diagrammatic representation is depicted on Fig. ||c. As has already been mentioned above, this term 
is generally spatially non-local. We take here only its local part. Technically it originates from a connected part 

of ^{S^2^\w,VK]si^\w,VK])^K, where iS*^^' 
cumbersome WW terms, one obtains 



is given by Eq. (116). Performing averaging over VK and omitting 



(z4^)z4^))vx = ^ ^.„..Tr{[M,T^^ + (M,T_J>^^_n(M-.,_^,,,rP„_n} 



7r^We^,e_ |^( 



(153) 





a) 



b) 




7^,(r, r) 



(154) 



In two dimensions this leads to the following logarithmic 
expression 



7^,(r,T) 



1 

Airg 



l/r 



duj 
2^ 



ln(r^'|c^|)F,_^(r,T), (155) 



c) 

FIG. o Diagrammatic represent ation o f the Gaussian part of 
the effective action, Sejf, Eq. (128b). a) and b) represent 
"out" and "in" relaxation terms c orre spondingly; they orig- 
inates from the first term in Eq. (149), {i S2)v K- Non-local 
term c) arises from the second term in Eq. (148), {iS2iS2)vK ■ 



For r2 = this expression coincide with the varia- 



tion of Bi^{r,T) over w in Eq. (145). Expanding to the 
first order in Q, one obtains the correction terms writ- 
ten in Eq. (pl3|). Eqs. (pl9|), (p^ and (p^) along 



with Eq. ( |113[ ) complete calculations of S'^iSeff /SwSw 
on the mean-field level. Let us notice for complete- 
ness that 5'^iSeff /SwSw = 0, which is a manifestation 
of the normalization condition. From another hand, 
d'^iSeff/SwSw 7^ 0, originates solely from the second 
term on the r.h.s. of Eq. ^^), cf. Fig. |c. 



B. Collisionless terms 

Finally we briefly discuss the physics of the collision- 
less terms. Collisionless terms originate from the real 
part of the selfenergy, TZ^^ (r, t) , and thus appear already 
in the first order in the bare interaction (unlike the colli- 
sion integral, which arises only in the second order). For 
the screened Coulomb interaction one obtains from Eq. 



(146) 



where we have used the superscript for the elastic mean 
free time, r*^', to avoid confusion with a physical time, 
T. If one linearize the kinetic equation ( |145| ) around the 
equilibrium distribution, its l.h.s. acquires the form 



Inf-i 
1 + — ^ 



T'='max{T, |e|} 



drF,{T)+ (156) 



d.F^i fde- 



A'Kg 



J^ln{r'^^\e-e'\)drFAr). 



We focus first on the logarithmic renormalization of the 
coefficient in front oidrF. This coefficient correspon ds to 
the charge Z in Finkel'stein's terminology |lj]. Eq. (15£) 
describes then the renormalization of Z (with the cor- 
rect coefficient). We stress, however, that in our theory 
renormalization of Z takes place at the level of the saddle 
point equation for the effective action, and not as a re- 
sult of the fluctuation corrections. This distinguishes Z 
from the conductance, g, whose renormalization occurs 
only at the level of the one loop correction, see section 
|VB| and Eq. (125). Physically renormalization of Z orig- 
inates from the suppression of the single-particle DOS 
by the residual short range interactions. This effect is 
due to the fact that single particle Hartree-Fock energies 
are shifted by the interactions in a way to reduce the 
DOS near the Fermi energy. We consider it very satis- 
factory that such purely mean-field effect is taken into 
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account by the saddle point equation and not by fluc- 
tuation corrections. The important thing, however, is to 
keep the last term of the expression ( |156| ) as well. This is 
to say that only the "out" -minus- "in" combination has 
the physical meaning. Being considered together, as an 
integral operator acting on i^^/, these terms do not lead 
to divergent corrections. 



VII. DISCUSSION 

We have developed a field theory for interacting dis- 
ordered metals using the Keldysh dynamic formulation. 
The advantages of this technique are twofold: (i) One 
avoids introduction of the replica trick; (ii) One natu- 
rally gains the ability to deal with non-equilibrium situ- 
ations. The latter manifests itself in the presence of the 
non-trivial object Ft^t'{r), which plays the role of the 
fcrmionic distribution function. The saddle point equa- 
tion of the theory turns out to be the kinetic equation 
which determines this function. No such object is ap- 
parent in the replicated Matsubara formulation Jl2| , p^ , 
since by construction it is limited to the equilibrium case. 
Based on the analogy with spin glasses |2^, one may 
speculate that non-trivial solutions Ft.t'{r) of the saddle 
point equation are analogous to the replica symmetry 
breaking solutions of the saddle point equation in the 
replica formulation. 

We mainly focussed our attention on the careful anal- 
ysis on the saddle equations of the theory. In particular 
we suggested the following two-step procedure: 

i) In the first step we account for the purely phase ef- 
fects of the fluctuating electric fields on the single particle 
Green function by an appropriately chosen gauge trans- 
formation. This enables us to get rid of the temporal 
variations of the Green function which are not related 
to the quasiparticle dynamics. The remaining temporal 
fluctuations of the Green functions are associated with 
the particle dynamics and can be described in terms of 
the quasiparticle distribution function Ft^t'{f)- This for- 
mulation ensures explicit gauge invariance of the kinetic 
equation and preserves the continuity relations at every 
stage. As a byproduct of this procedure we were able 
to obtain a non-perturbative expression for the DOS - 
the case where the phase effects give the main contribu- 
tion. Such phase effects do not contribute to gauge in- 
variant observables which are represented diagrammati- 
cally by closed loops. In the usual diagram technique this 
corresponds to a cancellation between certain diagrams 
(the diagrams containing double logarithms in two di- 
mensions) . By explicitly accounting for the phase effects 
we get rid of these diagrams which significantly reduces 
the number of terms in each order of perturbative expan- 
sion. 

ii) After the phase effects have been taken into account, 
we obtain a theory formulated in terms of the Q matrix 



field. The latter describes quantum fluctuations of the 
electron distribution function in the close vicinity of the 
Fermi energy. Restricting ourselves to the manifold of 
the massless fluctuations given by = 1 , we obtain the 
effective cr-model action, Seff[Q], Eq. (128b). Search- 



ing then for the extremum of this action, we arrive at 
the kinetic equation on the distribution function. After 
this two-step saddle point procedure one should consider 
the quantum fluctuations effects. The Altshuler-Aronov 
corrections to the conductance, g, turn out to be a 
manifestation of the one-loop quantum fluctuations. 

Although we have obtai ned r enorm aliza tion of both 
parameters g and Z, Eqs. ( 125 ) and (156), we deliber- 
ately avoided putting it in the framework of the renor- 
malization group. The point is that after introducing 
the phase transformation and integrating out the pho- 
ton fields, the effective action on Q, Eq. ( 128b| ), obtains 
a complicated form. We can not prove that this entire 
form is reproducible after the fast mode elimination. A 
seemingly better possibility is to perform the renormal- 
ization of the action, which contains both Q and V K 
fields, Eq. (^). In this case one has to specify how the 
relation between K and <I> fields changes in the process of 
renormalization. Since we believe that the introduction 
of the phase, K , is a vital element of the theory, the more 
complicated form of the action (compared to the one of 
Finkel'stein) is justifiable. 

We hope that the present formulation will help to shied 
light on the nature of the low-temperaturc phase of low- 
dimensional disordered metals. A few aspects of this the- 
ory seem to us very suggestive in this respect. Certain 
parallels with the spin glasses theory may prove to be 
useful. Apart from the extremely complicated problems 
relating to the character of the low-temperature phase, 
the functional Keldysh formalism may be useful for the 
description of non-equilibrium effects in disordered met- 
als and superconductors. The extension of this formalism 
to include the spin and Cooper channels will be a subject 
of our future work. 
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